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Abstract. In this article we consider viscous flow in the exterior of an 
obstacle satisfying the standard no-slip boundary condition at the surface 
, of the obstacle. We seek conditions under which solutions of the Navier- 

I Stokes system in the exterior domain converge to solutions of the Euler 

' system in the full space when both viscosity and the size of the obstacle 

' vanish. We prove that this convergence is true assuming two hypothesis: 

first, that the initial exterior domain velocity converges strongly in LP' to 
the full-space initial velocity and second, that the diameter of the obstacle is 
smaller than a suitable constant times viscosity, or, in other words, that the 
' obstacle is sufficiently small. The convergence holds as long as the solution 

' to the limit problem is known to exist and stays sufficiently smooth. This 

■ work complements the study of incompressible flow around small obstacles, 

' which has been carried out in [H El [6] 



oo 

o 



> ' 1. Introduction 

in 

\ The purpose of the present work is to study the asymptotic behavior of 

! famihes of solutions of the incompressible Navier-Stokes equations, in two 

^ ■ and three space dimensions, in the exterior of a single smooth obstacle, when 

both viscosity and the size of the obstacle become small. More precisely, 
let f2 be a smooth and bounded domain in M", n = 2, 3, such that Q is 
^ . connected and simply connected if n = 2 and M'^ \ f2 is connected and sim- 

■ ply connected if n = 3. Let e > and set = M" \ eQ. Let uq be a 

smooth, divergence- free vector field in M", which gives rise to a smooth so- 
lution u of the Euler equations, defined on an interval [0,T]. Let m'"'^ G 
L-((0,r);L2(n,)) n C'J[0,T);L\U,) n L^iiO,T);H^{Ue)) be a weak Leray 
solution of the incompressible Navier-Stokes equations, with viscosity z/, in 
He, satisfying the no-slip boundary condition at dUe- We prove that there 
exists a constant C = C{uo,fl,T) > such that if the following hypothesis 
holds: 



[H] We have that 



sup ||u'''^(-,0) - no||L2(n^^) 0, 

e<Cu 

as ^ 0, then we have that sup^^^,^ \\u'^''^ — it||Loo((o,T);L2(n5) ^ 0, as z/ ^ 0. 
Furthermore, if we assume that [H] occurs at a certain rate in z/ we can obtain 
an exphcit error estimate in L^. 

In addition, we prove that if we fix an initial vorticity uq in M^, smooth 
and compactly supported in \ {0} and consider u'^''^{-, 0) = A'^[u;o] + mif^, 
where denotes the Biot-Savart operator in Ilg, while is the normalized 
generator of the harmonic vector fields in n^, and m = J ujq, then hypothesis 
[H] is satisfied. In the case of the dimension three, if we fix an initial vorticity 
Uq in M^, smooth, divergence-free and compactly supported in M.^ \ {0} and 
consider u'^'^(-,0) = -^'^[ciJo]; where K"^ again denotes the Biot-Savart operator 
in n^, then it is proved in [1] that hypothesis [H] is satisfied. In both cases we 
have rates for the convergence of the initial data in such a way that Wu^'"^ — 
'^\\l^{Us) = C^(v^) when u 0, uniformly in time. 

A central theme in incompressible hydrodynamics is the vanishing viscosity 
limit, something naturally associated with the physical phenomena of turbu- 
lence and of boundary layers. In particular, a natural question to ask is whether 
the limiting flow associated with the limit of vanishing viscosity satisfies the 
incompressible Euler equations. This is known to be true in the absence of 
material boundaries, see [T71 [1] for the two dimensional case and [HI [20] for the 
three dimensional case. Also, if the boundary conditions are of Navier type, 
see [21 [161 (HI [23] , noncharacteristic, see [21] or for certain symmetric 2D flows, 
see [T8| [HI [15], convergence to an Euler solution remains valid. The most 
relevant case from the physical point of view corresponds to no slip boundary 
conditions. In this case, we have results on criteria for convergence to solu- 
tions of the Euler system, see [TUl [221 [231 [H] , but the general problem remains 
wide open. To be a bit more precise, let us assume that m° is a solution of 
the Euler equations in the exterior domain 11 and that u'^ is a solution of the 
Navier-Stokes equations with viscosity z/, with no-slip boundary condition in 
n. Suppose further that u'^ and m° have the same initial velocity Vq and that 

both M° and the family {u'^} are smooth, defined on a fixed time interval [0, T]. 
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It is easy to see that 

6E{u,t) = \\u%;t)-u''{-,t)\\LHTl) 

is uniformly bounded in u and t G [0,T] (by 2||t;o||L2 if ^^'o has finite energy), 
but it is not known whether 6E — > when u 0. In fact, given the experi- 
mentally and numerically observed behavior of high Reynolds number flows in 
the presence of boundaries, it is reasonable to conjecture that 6E does not, in 
general, vanish as ^ 0. Of course, this leaves open the possibility that u"^ 
might approach another solution of the Euler equations, different from u^. 

This article contains an answer to the following question: can we make SE 
small, by making both the viscosity and the obstacle small? This problem was 
one of the main motivations underlying the authors' research on incompressible 
flow around small obstacles. Our previous results include the small obstacle 
limit for the 2D inviscid equations, see [3, [13] and for the viscous equations, 
see [HI H]. The work we present here is a natural outgrowth of this research 
effort. 

The remainder of this article is divided into three sections. In Section 2, we 
state and prove our main result, namely the convergence in the small viscosity 
and small obstacle limit, assuming convergence of the initial data. In Section 
3 we study the problem of convergence of the initial data: for two space di- 
mensions, we adapt techniques developed in our previous work, while for three 
space dimensions we report to work by D. Iftimie and J. Kelliher. In Section 
4, we interpret the smallness condition on the obstacle as the condition that 
the local Reynolds number stays below a certain (small) constant. In addition, 
still in Section 4, we obtain an enstrophy estimate for the wake generated by 
the small obstacle and we list some open problems. 

2. Main theorem 

We use the notation from the introduction to state and prove our main 
result. We consider the initial- value problem for the Navier-Stokes equations 
in Ilg, with no-slip boundary condition, given by: 

' dtu"'' + u""'' ■ Vu"^' = -Vp"'' + ^/An''■^ in Yi, x (0, oo) 

div u"'^ = in X [0, oo) 

I M'''^(x,t) = forxeSn^, t>0 

^ u''^^{t = 0) = §^{x) for X e n^, t = 0. 
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We assume that the initial velocity 'd^ G L'^^JJlg) is divergence-free and 
tangent to dlle-, but we do not assume that it satisfies the no-slip boundary 
condition. In three dimensions we assume further that G L'^ijle). Under 
these hypothesis' it was shown by H. Kozono and M. Yamazaki, see [I2], that, 
in two dimensions, there is a unique global strong solution to ([T]) with initial 
velocity d^, while, in three dimensions, there is a global Leray weak solution 
u^''^ of ([T]), see [3]. More precisely, in three dimensions there exists 

M^'^ G L°°([0, oo); L2(n,)) n C° ([0, oo); ^^(n,)) n L2^,([0, ^)-Hl{Yi,)) 

such that u^'^ is a distributional solution of ([1]) and the following energy in- 
equality holds true: 

(2) ||n'^'^(t)|U^(n,) + / || Vn'^'^(s)|U2(n,) ds < H\\L^(n.) Vt > 0. 

Jo 

Both and are defined only in He, but we will consider them as 

defined on the whole space by extending them to be identically zero inside eQ. 

Let Mo be a smooth, divergence-free vector field defined in all M", and let u be 
the corresponding smooth solution of the Euler equations; in two dimensions 
u is globally defined while in three dimensions it is defined, at least, on an 
interval [0, T]. 

We are now ready to state our main result. 

Theorem 1. Assume that 

W - %IU2(Mn) ^ as e ^ 0. 

Fix T > 0, arbitrary if n = 2, and smaller than the time of existence of the 
smooth Euler solution ifn = 3. Then there exists a constant Ci = Ci{Q, Uq, T) 
such that, if e < Civ, then 

- n(-,t)||i2(iKn) as ^ 0. 

Moreover, if we assume that H"*?^ — mo||l2(]r.i) = 0{^), then there exists also 
C2 = C2{T, uq, Q) such that 

for alio <e < Ciu and allO<t<T. 
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Before we proceed with the proof, we require two technical lemmas. To state 
the first lemma we must introduce some notation. As in the statement of the 
theorem, u denotes the smooth Euler solution in M."^. In dimension two, we 
denote hj ip = il){x,t) the stream function for the velocity field u, chosen so 
that ^^(O,^) = 0. In dimension three, ip denotes the unique divergence free 
vector field vanishing for x = and whose curl is u. In other words, we set 

in dimension two so that u = V'^ip and 

in dimension three so that u = curl-?/;. In both two and three dimensions one 
has that ip and Vip are uniformly bounded on the time interval [0,T]. 

Let i? > be such that the ball of radius R, centered at the origin, contains 
Q. Let if = (p{r) be a smooth function on M.,. such that ip{r) = if < r < 
R+1, ip>0 and <^{r) = 1 if r > i? + 2. Set = <^'{x) = <^{\x\/e) and 

in dimension two and 

in dimension three. In both dimensions two and three, the vector field is 
divergence free and vanishes in a neighborhood of the boundary. 

We also re-define the pressure p = p{x, t) from the Euler equation in 
with data uq so that p(0, t) = 0. 

Lemma 2. Fix T > 0. There exist constants Ki > 0, i = 1, . . . ,5 such that, 
for any < e < Eq and any < t < T we have: 



(1) 
(2) 
(3) 
(4) 
(5) 



Im^II^oo < K2, 

|V^V<^^||loo + IIV^VVIU- < K,/€, 
\pVip'\\L2 + \\dtiljVip'\\L2 <K^e. 



Above, we used the notation || V^/'Vv2'^||l°o = 'Yliij \\'^i'4^djip^\\L°° in dimen- 
sion two and || V-j/^Vv^^Hloo = Ylijk \\9i'ipkdjf'^\\L°° in dimension three. Similar 

notations were used for the other terms. 
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Proof. Some of the inequalities above can be improved in dimension three. 
However, it turns out that these improvements have no effect on the statement 
of Theorem [H Therefore, to avoid giving separate proofs in dimension three 
we chose to state these weaker estimates. 

Recall that both u and Vu are uniformly bounded. First we write 

in dimension two and 

= di cnT\{ip'^'ip) = udicp'^ + Vv9^ x diip + diV^p^ x t/) + ip'^diu 

in dimension three. The supports of the first three terms of the right-hand 
sides of the relations above are contained in the annulus e{R + 1) < |x| < 
e(-R + 2), whose Lebesgue measure is C(e^). Furthermore, |Vv?^| = 0{l/e), 
IVV^I = 0{l/e'^) and \il){x,t)\ = 0{e) for |x| < e{R + 2), since ip{0,t) = 0. 
Taking norms in the expressions above gives the first estimate. 

Next we observe that = (f'^u + ipV'^ip'^ or = ip'^u + Vy^^ x ip. Clearly 
(f'^u is bounded and to bound the second term, we use again that ip{0,t) = 0, 
which proves the second estimate. For the third estimate, observe that u'^ — u 
and — ip^u are bounded, as we have just proved, and have support in the 
ball |x| < e{R + 2). For the fourth estimate, we use again that '?/'(0,t) = 0. 
The last estimate follows from two facts: that the functions whose L^-norm 
we are estimating have support on the ball |x| < e{R + 2) and that they are 
both bounded, since p(0, t) = and V'tlO, t) = 0. □ 

We also require a modified Poincare inequality, stated below. This is fairly 
standard, but we include a sketch of the proof for completeness. 

Lemma 3. Let 1^ he the obstacle under consideration and let R he such that 
Vt C Bfi. Consider the scaled obstacles eVt and the exterior domains II^. Then, 
ifWe H^iUe) we have 



L2(n,nB(fl+2)j < KQe\\VW\\L2(u,nB^a+2)e)- 

Proof. The proof proceeds in two steps. First we establish the result in the case 
e = 1. Suppose, by contradiction, that there exists a sequence {W^'^} C Hq(IIi) 
such that \\W''\\L2(u,nBii+2) > k\\'^W'\\L2(u^nBn+2)- Set 



\W'\\L2(mnBn+2)' 



Then G ifg (Hi), with unit norm in Hi fl while the L^-norm of its 

gradient on Hi fl -B_r+2 vanishes as k oo. Thus, passing to a subsequence if 
necessary, V'' V, weakly in and strongly in (on Hi fl Bf{^2) so that 
||^||L2(ninBH+2) ~ ^^'^ ~ 0- Since this set is connected, it follows that V 
is constant in Hi ni?/j+2- By continuity of the trace, the trace of on F = dUi 
must vanish, which shows that V = in Hi fl -Bi?+2, a contradiction. 

We conclude the proof with a scaling argument. Let W G H^(Il^) and set 
Y = Y{x) = W{ex). Then Y G Hq{IIi). Using the first step we deduce that 
there exists a constant Kq such that ||^||L2(ninBi{+2) — ^6\\^'^\\L'^iUinBn+2)- 
Undoing the scaling we find: 

\\wF 

\W{ex)\'^dx = - 



\Y\\Un.nB,,2, = I mex)\^dx = " 



l|Vr|li.(n,ni^.,2) = / e'\VW {ex)\' dx = e'-'^WVWWl.^^^^^ y 

JUinBa+2 

The desired result follows immediately. 

□ 

We are now ready to prove Theorem [H 

Proof of Theorem [Tt We begin by noting that, since u is a smooth solution 
of the Euler equations in M" x [0,T], it follows that 

(3) = H-,t)\\L2^,n)<Ce<CV^, 

ii e < Cu. Hence it remains only to estimate the L^-norm of the difference in 
Ilg, which we do below. 

We first give the proof in two dimensions and then we indicate how the proof 
should be adapted to three dimensions. 

2.1. Case n = 2. The vector field u'^ is divergence free and satisfies the equa- 
tion 

We set Ty^'^ = u"''^ — . The vector field W^'^ is divergence free, vanishes on 
the boundary and satisfies: 
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We perform an energy estimate, multiplying this equation by W"'^ and in- 
tegrating over He. We obtain 



(4) —WW-^lh + uiivw^^^wl, 



'He 



[ ■ [(m"^'^ ■ V)u''^'] dx+ [ W^' ■ [{^'u ■ V)u] dx 

+ / -^'Vpdx- I W^' ■ dti'V^ip' dx. 



We will examine each one of the five terms on the right-hand-side of iden- 
tity (jll). We look at the first term. We use Cauchy-Schwarz and Young's 
inequalities followed by Lemma [21 item (1), to obtain 



(5) 



u I VW'' ■ Vu' dx 



< -( \\VW''^'\\l2 + Ki ]. 



Next we look at the second and third terms together. We write 

- / ■ [{u"^' ■ V)u'^'"] dx+ I ■ [((^% ■ V)n] dx 

- I ■ [{{W^' + u') ■ V)(iy"'" + u')] dx+ I ■ [{^'u ■ V)n] dx 

- I ly"'" ■ [(ly"'" ■ v)u'] dx- I ly"'" ■ [{u^ ■ v)u'] dx 

+ [ W'' ■ [{ip'u ■ V)n] dx , 

where we used the fact that J W'" ■ [((VT'^'^ + m^) ■ V)^'^] = 0. Finally, we 
add and subtract / W'^''^ ■ [{u'^ ■ V)u] to obtain 



IXI 



< 



[ W'' ■ {W^' ■ Vm") dx + [ W^' ■ [u' ■ V(m - u')] dx 

+ / W"^' ■ [{{^'u - u') ■ V)u] dx 
I -{W^' ■Vu')dx + [ {u-u') -[{u' ■\/)W'''']dx 

I W"^' ■ [{{^'u - u') ■ V)u] dx 



< 



/ {W^' ■Vu')-W''^'dx +\\u-u'\\l2\\u'\\lo^\\VW''^'\\l2 

+ \\W''''\\^2y-U-U%LA\^u\\Lo 



For each z = 1 , 2 we have that 

(6) diu' = di^V^if' + ipdiV^ip' + diip'u + ip'diU. 
Therefore, 

(7) |j| < (||v^AV<^^|u^ + im^^mL^^W'ml,^^^^ + y'Vu\\L^\\w''''\\l, 

+ ||n - u'\\l2\\u'\\loo\\VW''^'\\l2 + \\W''\\L2yu - ^^1^2 1| Vn||i.o, 

where is the set fl i?(R+2)e, which contains the support of Vip'^. 

Hence, using Lemma [2], items (2), (3) and (4), together with Lemma [3l in 
the inequahty ([Tj), we find 



(8) |X| < ^-je'Kl\\VW'mi2 + y'^uWLooWW'mi^ 

+ K2K3e\\VW''\\l2 + K3e\\Vu\\Loo\\W'mL2. 

Next we look at the fourth and fifth terms in (jlj). Recall that we chose the 
pressure p in such a way that p(0, t) = 0. We find 



\J\ 



< 



W'' ■ pWif' dx 



■ dtipV^^' dx 



We estimate each term above to obtain, using Lemma [2] item (5), 
(9) \J\ < {\\pVv'\\l2 + ll^t^VVlUOIlW^'^'IlL^ < i^5^||Vr^'1|L2. 
We use estimates ([H]) and © in the energy identity (jl]) to deduce that 



1 d 



+Ko\\W 



L2 + 



Above we have used the notation Kq = sup^ Hv^'^VuHlcx) and ^5 = ^5 + 
-f^sll Vm||loo. 

At this point we choose e so that 

(11) 

With this choice, letting y = y{t) = ||Vr'^'^||^2, we obtain 

(12) § < + ^KiKf-^ + Ky + {2K, + l)y < C[v + C',y. 
From Gronwall's inequahty it follows that 

(13) Wu"'' - u^Wl^^^,) < CiT,uo,n) (/^+ ||r-ng||i.(n.)) • 

If ll??^— uo||l2 then it follows from (fT3l) together with LemmaEl item (3), 
and from ([3]), that \\u'^''^ — u\\l2 — > 0, as desired, where the constant Ci can be 
chosen to be {8K4^Kq)~^. If we assume further that H^?^ — mo||l2 = C^(v^) then 
the second part of the statement of Theorem [T] easily follows. This concludes 
the proof in the two dimensional case. 

2.2. Case n = 3. The proof in dimension three is similar to the previous one. 
There are two differences: notation and the justification that we can multiply 
the equation of W'^ by W'^. 

First, about notation. One has to replace everywhere the term dtipV^ip^ by 
V(/?^ X dtip and also relation becomes 

These two modifications are just changes of notations. These new terms are 

of the same type as the old ones, so the estimates that follow are not affected. 

Second, we multiplied the equation of W"'^ by W^'^. The solution u"'"^, and 

therefore W'' too, is not better than L°°(0, T; ^2(0^)0^^(0, T; H\Ue)). But 

it is well-known that some of the trilinear terms that appear when multiplying 

the equation of W'^''^ by W^'^ are not well defined in dimension three with this 

regularity only. In other words, one cannot multiply directly the equation of 

W^'^ by W^'^. Nevertheless, there is a classical trick that allows us to perform 

this multiplication if the weak solution u'^'^ verifies the energy inequality. What 

we are trying to do, is to subtract the equation of from the equation of n'''^ 

and to multiply the result by m'''^ — m^. This is the same as multiplying the 

equation of m'''^ by u'^''^, adding the equation of times and subtracting 
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the equation of m^'^ times u'^ and the equation of times m^'^. Since is 
smooth, all these operations are legitimate except for the multiplication of the 
equation of u'^'^ by u'^'^. Formally, multiplying the equation of u^''^ by n'''^ 
and integrating in space and time from to t yields the energy equality, i.e. 
relation ([2D where the sign < is replaced by =. Since we assumed that the 
energy inequality holds true, the above operations are justified provided that 
the relation we get at the end is an inequality instead of an equality. But an 
inequality is, of course, sufficient for our purpose. Finally, to be completely 
rigorous, one has to integrate in time from the begining. That is, we would 
obtain at the end relation ( fT2|) integrated in time. Clearly, the result of the 
application of the Gronwall lemma in (fT2l) is the same as in (fT2l) integrated in 
time. This completes the proof in dimension three. 

Remark: The proof above is closely related to the proof of Kato's criterion 
for the vanishing viscosity limit in bounded domains, see [TU]. Both results are 
based on estimating the difference between Navier-Stokes solutions and Euler 
solutions by means of energy methods. In Kato's argument, the difference is 
estimated in terms of the Navier-Stokes solution, on which Kato's criterion 
was imposed. In contrast, our proof estimates the difference in terms of the 
full-space Euler solution, which is smooth in the context of interest. 

3. Compactly supported initial vorticity 

Now that we are in possession of Theorem [T] we will examine two asymp- 
totic problems for which we can prove the convergence condition on the initial 
velocity. We focus on flows with compactly supported vorticity, and the diam- 
eter of the support of vorticity becomes the order one length scale, relative to 
which the obstacle is small. 

Let us begin with the three dimensional case. We consider an initial vorticity 

ujq which is assumed to be smooth, compactly supported in M? \ {0}, and 

divergence-free. Let e > be sufficiently small so that the support of ujq is 

contained in He. The domain He is assumed simply connected so that there 

exists a unique divergence-free vector field, tangent to c^II^, in L^(ne), whose 

curl is cjq, see, for example, [1]. We take d'^ to be this unique vector field. We 

take Mo to be the unique divergence-free vector field in L^(M^) whose curl is 

cuo, given by the full space Biot-Savart law. 
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In D. Iftimie and J. Kelliher studied the small obstacle asymptotics 
for viscous flow in U^, for fixed viscosity, in three dimensions. They proved 
that the small obstacle limit converges to the appropriate Leray solution of 
the Navier-Stokes equations in the full space. One important ingredient in 
their proof was precisely to verify strong convergence of the initial data; in our 
notation, Iftimie and Kelliher proved that 

- %IU2(K3) = 0(et), 

as £ — > 0. 

We may hence apply Theorem [1] to obtain the following corollary. 

Corollary 4. Let uq G C^(M'^\{0}) and consider uq and'd'^ defined as above. 
Fix T > and assume that the solution u = u{x,t) of the incompressible 
Euler equations in , with initial velocity Uq, exists up to time T. Let u'"'^ 
be a Leray solution of ([l]) with initial velocity . Then there exist constants 
Ci = Ci{Q, Uq, T) > and C2 = 6*2(^2, ci^o, T) > such that 

for allO<e < Civ and allO<t<T. 

Next we discuss at length the case n = 2. In dimension two the exterior do- 
main is no longer simply connected. This means that the vorticity formulation 
of the Euler equations is incomplete, and we must specify the harmonic part 
of the initial velocity as well as the initial vorticity, see [H] for a thorough dis- 
cussion of this issue. To specify the asymptotic problem we wish to consider, 
we must choose the initial data for ([T]). 

Let Uq be smooth and compactly supported in \ {0}. Let K denote 

the operator associated with the Biot-Savart law in the full plane and set 

H = x^/(27r|xp), to be its kernel. Let K'^ be the operator associated with the 

Biot-Savart in 11^, i.e., = V'^A^^, where An^ is the Dirichlet Laplacian in 

Ilg. Let LL^ be the generator of the harmonic vector fields in 11^, normalized 

so that its circulation around dUs is one. The divergence-free vector fields in 

with curl equal to uq are of the form /^'^[cuo] + oH'^, with a G M, see [S]. 

In [6] the authors studied the asymptotic behavior, as e ^ 0, of solutions of 

([T]) with u fixed and initial velocity /'^'^[co'o] + It was shown in [U] that 

m'''^ converges to a solution of the Navier-Stokes equations in the full plane 
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with initial data K[ujq] + (a — m)H, where m = J ujq, as long as |a — m| is 
sufficiently small. 

For the vanishing viscosity limit, we must consider only the case a = m. 
There are two reasons for this. First, K'^[ujQ]+aH'^ converges weakly to -ft'[co'o] + 
(a — m)iir in distributions, see Lemma 10 in [Hj, but, as we shall see, this 
convergence is not strong in (see Remark 1 following the proof of the next 
lemma). Second, one cannot expect solutions of the Euler equations in the full 
plane with initial velocity Kftuo] + (« — m)H to be smooth (even existence is 
not clear) unless a = m. 

In view of this discussion, set uq = K[ujq] and = /'^^[ciJo] + mH^. With 
this notation, we can prove strong convergence in of the initial data, as 
follows. 

Lemma 5. Fix Eq such that the support of ujq does not intersect fi^ for any 
e < Sq. There exists a constant C > 0, depending on Vt and ujq such that 

W - ^o|U2(M2) < Ce. 

Proof. We begin the proof with a construction whose details can be found in 
[5]. In Section 2 of [5], an explicit formula for both K'^ and H'^ can be found 
in terms of a conformal map T, which takes 11 into the exterior of the unit 
disk centered at zero. The construction of T and its behavior near infinity are 
contained in Lemma 2.1 of [5]. Using identities (3.5) and (3.6) in [5], we have 
that the vector field H'^ can be written explictly as 

and the operator can be written as an integral operator with kernel /C^, 
given by 

^e^^r^ru ,.( inx/e)-T{y/e))^ {T{x/e) - {Tjy/e))*)^ 
2ne ^ ^ '\ \T{x/e)-T{y/eW \T{x/e) - {T{y/e))*\^ 
where x* = denotes the inversion with respect to the unit circle. Fur- 

thermore, we recall Theorem 4.1 of [5], from which we obtain 

Wi^^h^iu,) < C\\uJo\\]^^\\uJo\\]!i, 

for some constant C > 0. 

To understand the behavior for e small in the expressions above, we need to 
understand the behavior of T{x) for large x. We use Lemma 1 in [7j, which is a 
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more detailed version of Lemma 2.1 in [5], to find that there exists a constant 
P > such that 

(14) T{x/e) = /3xe~^ + h{x/e), 

with h = h{x) a bounded, holomorphic function on IIi satisfying \Dh{x)\ < 
C/|xp. Therefore, 

(15) \DT{x/e)- (31\<C- 



w - 


- Z2\ 




N2I 



We will need a further estimate on the bounded holomorphic function h 
h{z) = T{z) — Pz, namely that 

(16) IHz,) - h{z2)\ < C 

for some constant C > independent of Zi, Zi- This estimate holds since, 
by construction (see Lemma 2.1 in [5]), we have that hi^z^ = g{l/z) with g a 
holomorphic function on (Hi)*, whose derivatives are bounded in the closure 
of (Hi)*. Here, (Hi)* denotes the image of Hi through the mapping x ^— x* = 
x/|a;p to which we add {0}. Here, we are using the following fact: If D is 
a bounded domain with dD a Jordan curve then any bounded function 
/ G C^{D) with bounded derivatives is globally Lipschitz in D. This fact is a 
nice exercise in basic analysis, which we leave to the reader. 
Therefore we have 



(17) \h{z^) - h{z2)\ 







< c 


1 


1 


















Z2 



c 



\zi - 


- Z2\ 


\zi 


\Z2\ 



In order to estimate \\'d^ — Uq\\l2(ji^^ we use the fact that the support of ujq 
is contained in H^ for e sufficiently small e to write 

^-i^''-)-"°wi-l(;^^-'-/-' ^ri:;tri:;? -i^)-°'^)^^ 

-e"^ <"/^>(MVi)F- |r(./.) - (r(,/.)).|= h''^^'^ = « ■ 

Let us begin by estimating *B^. We make the change of variables 77 = 
eT{y/e), whose Jacobian is J = | det{DT~^) {r] / e)\, a bounded function. Ad- 
ditionally, we set z = eT{x/e). With this we find: 

W = DT\x/e) [ - ^-^-^^UieT-'iv/e)) Jdr^. 



14 



We observe now that there exists p independently of e such that the support 
of ujQ{eT~'^{ri/ e)) is contained in the set {|?7| > p}. Moreover, one can bound 
\z — e'^rj*\ > \z\—e'^\ri*\ > \z\ —e^/p > \z\/2 provided that < p/2. Therefore 
we can write 

.2 



\w\ < c 



^'1 


1 * 1 




\z\ 


\z - 





e 



\z\ 

where C depends on the support of cuq, on the L^-norm of ujq and on the 
domain VL through the bounds on the conformal map T and its derivatives. 
Finally, we use this estimate in the integral of the square of QS*^: 

dx <Ce^ I -X7dz< Ce^, 
J{\z\>e} \z\ 

as desired. 

Next we treat 21^. First we re- write 21^ in a more convenient form: 

" = y„_ (-/-)l7 |r(V.) - T(yM\-' - ^I^Hy^'« 

By ( ITSll . the term 21^2 can be easily estimated: 

|2t^2|<^ / 1 :\uo{y)\dy<C 



1^ ~ y\ ' 

so this reduces to an estimate similar to the one we found for 53^. 

Next we examine 2t^i. We use the expression for T given in fll4l) to write 

;f)Kf)-Mf; 



With this we have: 



\^\hi^)-hm]\ 



In. l^-2/ + (f)[Mf)-MPII^-y| 

We will make use several times of the estimate we obtained for h given in (fT71) . 



First 
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Using (fT8|) gives 

Let -R, r > be such that the support of ujq is contained in the disk of radius R 
and outside the disk of radius r. We will estimate in two regions: |x| > 2R 
and |x| < 2R. Also, recall that the estimate of the L^-norm of is to be 
performed in 11^ so we may assume throughout that |x| > Ce. Suppose first 
that |x| > 2R. Then we find: 

Wi\<c—. 

Above we used that r < \y\ < |x|/2 and hence \x—y+{^)[h{^) — h{'^)]\ > C|x| if 
e is sufficiently small, since h is bounded. Finally, in the region Ce < |x| < 2R 
we use ( !T6l) and the fact that \y\ is of order 1 to bound 

\^-y + ) - Mf )] > ll^; - y| - e\\x - y\/\x\\y\)\ > 

for e small enough. Therefore 

\^M<C^ ! P^dy<C^<Ce. 

fI F ^ y\ fI 

Clearly this last portion has L^-norm in the region |x| < 2R bounded by Ce. 

□ 

Remark 1: Let ujq G C^{R^ \ {0}) and a G M, a 7^ m. We observe 
that A'"[ti;o] + (^H^ does not converge strongly in to -R'fct'o] + (a — m)H. 
We argue by contradiction, assuming this convergence holds. In view of the 
Lemma above, i^"[u;o] + mH^ converges strongly to -^'[(^o] so we must have 

{a — m)H' (a — m)H. 

This does not hold, as it can be easily seen in the case of the exterior of the 
disk. In this case, if" = H outside the disk of radius e, but vanishes for 
|x| < e. Since ||if||L2({|x|<£}) = 00, we have a contradiction. 

Remark 2: Note that if we were willing to confine our analysis to the 
exterior of a small disk, the proof of Lemma [5] would be much simpler. Indeed, 
let fie = {|x| > e}. Then the conformal map T is the identity, so 21" = and 
all that is needed is the easier estimate for 03". 
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Remark 3: The constant a — m is precisely the circulation of ■i?'^ around 
the boundary of 11^. The condition a — m = is physically reasonable, in 
particular because viscous flows vanish at the boundary, and therefore, so 
does their circulation. This is the condition for the small obstacle limit of 
ideal flow to satisfy Euler equations in the full plane, see [3] and also for the 
small obstacle limit of viscous flows to satisfy the full plane Navier-Stokes 
equations for all viscosities, see [0]. The argument in [H] required sufficiently 
small a — m to obtain the appropriate limit when e —>■ 0, and the smallness 
condition was actually a — m = 0{v) as — > 0. 

We conclude this section with the formal statement of a corollary which 
encompasses Theorem [T] and Lemma [5l 

Corollary 6. Let ujq G C^iM? \ {0}) and consider uq and defined as in 
Lemma\^ Let u = u{x,t) be the global smooth solution of the incompressible 
Euler equations in M^, with initial velocity Uq. Let u^'^ be the solution of ([T]) 
with initial velocity . Fix T > 0. There exist constants Ci = c^o, T) > 

and C2 = C2{fl,uJo,T) > such that 

for allO<e < Civ and allO<t<T. 

4. Physical interpretation and conclusions 

The behavior of incompressible viscous flow past a bluff body, such as a long 
cylinder or a sphere is a classical problem in fluid dynamics, to the extent of 



having conference series devoted to it, see http://www. mae. comell. edu/bbviv5/. 
Let us consider the simplest situation, two-dimensional flow of a viscous fluid 
with kinematic viscosity /i, flUing the whole plane minus a disk of diameter L, 
with constant driving velocity U at inflnity. The disturbance caused by the 
disk, known as its wake, depends only on the Reynolds number associated with 
the flow, given by 

/i 

The observed behavior of the wake begins, for small Re, as a steady solution of 
the Navier-Stokes equations, but the wake undergoes a series of bifurcations as 
Re grows, progressively developing steady recirculation zones (4 < i?e < 40, 
periodic recirculation and a Von Karman street (40 < Re < 200), nonperiodic 
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vortex shedding (200 < Re < 400), leading to turbulence {Re > 400). See 
[To] . Section 15.6, for details and illustrations. 

In our problem, which involves nearly inviscid flow past a small bluff body, 
the qualitative behavior of the wake of the small obstacle is determined by the 
local Reynolds number, which encodes the way in which an observer at the scale 
of the obstacle experiences the flow. Basically, by making our obstacle small, 
we are making the flow more viscous at its scale. We assume that the Navier- 
Stokes system under consideration, ([T]), is dimensional, i.e. has time and space 
measured in seconds and meters, and mass normalized so that fluid density is 
one. In these units, the kinematic viscosity for air is // = 14.5 x 10~^m^/sec, 
and for water it is /i = 1.138 x 10~^m^/sec, both at 15°C. 

Let us restrict our discussion to the two-dimensional case. The smallness 
condition in Theorem [1], (fTTl) . reads e < Cifi, and the dimensional constant 
Ci, requires closer scrutiny. Actually, the constant Ci is given by: 

1 



where K4 appears in Lemma [21 item (4), and Kq is from Lemma El Kq is a 
non-dimensional constant that depends on the shape of the obstacle Q. The 
constant K4 can be chosen as 

K^= sup e{\V^lJVlp'\{x,t) + \xljV^ip'\{x,t)). 

The function tp above is the stream function of the full-plane Euler flow, 
adjusted so that ip{0,t) = 0. Also, eVip'^ is C(l), localized near the obstacle 
and eV'^ip'^ is 0{l/e), also localized near the obstacle. Therefore, both terms 
included in K4 are associated with flrst derivatives of the stream function 
at the obstacle, i.e. with the local velocity u{0,t). Therefore, we can write 
K4 = supj |u(0, t)| (we can also assume that the limiting Euler flow is 
stationary, to avoid the time dependence). 

From the point of view of the obstacle, the inviscid velocity u{0,t) acts as 
a constant (in space) forcing velocity imposed at inflnity, and therefore, the 
qualitative behavior of the wake of the obstacle is determined by the local 
Reynolds number Rei^^ = u{0,t)e/fi. Clearly, condition (fTTj) can be rewritten 
as 

1 



ReiQc < 



8K^K^ 
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The non-dimensional constant is related to the constant in the Poincare 
inequality in the unit disk. Examining our proof for the case of the disk, we 
cannot make the constant K^Kq smaller than something of the order of 10. 
Therefore, our result is restricted to rather viscous wakes. 

When it occurs, the turbulence is caused by vorticity shed by the obstacle 
through boundary layer separation. The main difficulty in studying the van- 
ishing viscosity limit in the presence of boundaries is the fact that, although 
the Navier-Stokes equations do have a vorticity form, valid in the bulk of the 
fluid, the vorticity equation does not satisfy a useful boundary condition, so 
that we cannot control the amount of vorticity added to the flow by the bound- 
ary layer. In the proof of Theorem [T], we found a way of controlling the kinetic 
energy of the wake without making explicit reference to the vorticity. At this 
point, it is reasonable to ask whether we can control the vorticity content of 
the wake as well. To answer that, we introduce the enstrophy of the 

flow: 



Of course, enstrophy measures how much vorticity is in the flow, but its be- 
havior as z/ — s> is also involved in the statistical structure of a turbulent 
wake. 

Corollary 7. For any T > there exists a constant C > 0, independent of v 
such that 



Jo 

Proof. We go back to relations (ITOl) and ( fT2l) and include the viscosity term 
which had been ignored. We flnd: 





+ ^\\VW'''%,<C[u + C',y. 



We next integrate in time to obtain 




^0 

Now we use Theorem [1] and ignore a term with good sign to obtain 
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where we used Lemma |5] together with item (3) from Lemma [2] to estimate the 
initial data term. From this we conclude that 

/ \\\/W''^'{-,t)\\l2dt<C. 
Jo 

Finally, we observe that 

n"'' < cwvw^ml^ + c\\Vu'\\lo_ < cwviV'^wl, + ck,, 

by item (1) in Lemma [21 This concludes the proof. 

□ 

Finally, let us consider some open questions naturally associated with the re- 
search presented here. First, one would like to weaken, and ultimately remove, 
the smallness condition on the size of the obstacle; this is the most physically 
interesting follow-up problem. Second, one would also like to consider two 
dimensional flows with nonzero initial circulation at the obstacle, in order to 
study the interaction of the vanishing viscosity and vanishing obstacle limits 
in more detail. This would improve the connection of the present work with 
the authors' previous results in [HI E] • An easier version of this second prob- 
lem would be to consider an initial circulation of the form 7 = 7(2/) and find 
out how fast 7 has to vanish as z/ in order to retain our result. A third 
problem is to describe more precisely the asymptotic structure of the difference 
between the full-space Euler flow and the approximating small viscosity, small 
obstacle flows. 
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